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Abstract 

It is shown that the observed anisotropy of magnetization at high magnetic 
fields in RbMnBra, a quasi-one-dimensional antiferromagnet on a distorted 
stacked triangular lattice, is due to quantum and thermal fluctuations. These 
fluctuations are taken into account in the framework of linear spin-wave the- 
ory in the region of strong magnetic fields. In this region the divergent one- 
dimensional integrals are cut off by magnetic field and the bare easy-plane 
anisotropy. Logarithmical dependence on the cutoff leads to the "enhance- 
ment" of the anisotropy in magnetization. Comparison between magnetiza- 
tion data and our theory with parameters obtained from neutron scattering 
experiments has been done. 

I. INTRODUCTION 



The question of the effect of fluctuations on the magnetic properties of quasi-one- 
dimensional antiferromagnets has been intensively discussed in the last few years. Particular 
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attention has been given to the cases where there is a stacked triangular lattice and an an- 
tiferromagnetic interaction, so that there is a frustration in going from one-dimensional to 
three-dimensional ordering. 

This situation is apparent in materials with hexagonal CsNiCls-type crystal structure 
having the general formula ABX3 where A is an alkali metal, B is a bivalent metal of 
3d group, and X is a halogen. While spin dynamics of CsNiCla and RbNiCls (S=l, easy- 
axis anisotropy) has been studied to examine Haldane conjecture the compounds of 
the same group, CsMnBrs, RbMnBra (S=5/2, easy-plane anisotropy) can be considered as 
appropriate for checking the applicability of standard spin-wave theory for the case of large 
half-integer spins. 

CsMnBrs has been investigated by various experimental techniques as an example of a 
frustrated antiferromagnet on a stacked triangular lattice exhibiting a number of unusual 
magnetic properties, such as a field-induced phase transition from the triangular phase 
to the collinear one, critical behavior associated with a chiral degeneracy and an unusual 
phase diagram. In addition to numerous neutron scattering experiments and to ESR 
measurements [Q, the magnetization process has been studied in detail with results 
which are in good agreement with classical calculations ||^ except for two discrepancies. 
Namely: 

• the measured magnetic torques are significantly smaller than the theoretical ones; 

• in large magnetic fields {H > He, where He is the critical field of the transition to 
the collinear phase) there is a considerable anisotropy between magnetization when 
magnetic field is applied along C-axis of the crystal (M^"*") and the one (M^-*"*") for 
magnetic field along easy plane absent in the theory 

It was suggested that in order to describe these peculiarities the anisotropy of quantum 
fluctuations should be taken into account. In this paper we consider on the basis of another 
compound RbMnBrs how both quantum and thermal fluctuations affect the magnetization 
of quasi-one- dimensional noncoUinear antiferromagnet. 
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II. CRYSTAL STRUCTURE AND MAGNETIC PROPERTIES OF RbMnBrg 



Powder neutron scattering measurements by Glinka et al. M showed that RbMnBra 



exhibits antiferromagnetic order below = 8.8±0.1K. Single crystal measurements |10.11 



show that the magnetic structure of RbMnBrs is incommensurate with Mn"*"^ moments lying 
in the basal plane and with antiferromagnetic ordering in the c direction. The commensurate 
ordering in CsMnBra gives scattering peaks of type (|, |, 1) [0] while in the incommensurate 
RbMnBrs these peaks are replaced by two triads of peaks near to (|, |, 1) with three fold 
symmetry around (|, |, 1). Because the incommensurate peaks surround a commensurate 
position in reciprocal space, it is likely that the incommensurate structure is locally similar 
to the 120° triangular structure of CsMnBrs, but with angles between neighboring spins in 
the basal plane somewhat larger than 120°. 

The appearance of incommensurate magnetic structure may be due to the presence of 



crystal structure distortions. According to X-ray study by Von Fink and Seifert [0] as the 
temperature decreases the crystal structure of RbMnBrs changes at 470K from CsNiCls- 
type structure (space group PG^/mmc) to the so-called KNiCls-type structure (space group 



PQ^cm). Measurements of the birefringence [|T3| and recent neutron scattering experiments 



1^ show that another structural phase transition takes place at about 220K. As reported 
13 below this temperature the crystal has orthorhombic structure with cell dimensions 



m 



a=14.680A, 6=12.805A, c=6.516A at 12K. 

At low temperatures a field of 3.0T applied along the a axis produces a first-order phase 
transition to a commensurate structure |[TT1JT5| , so that the lattice parameter of the magnetic 



cell is eight times that of the nuclear cell ||I0[. In a magnetic field of about H^. = 3.9T and 
above applied in the basal plane, the magnetic structure becomes collinear, as in the case 
of CsMnBrs with He = 6.4T. This conclusion can be derived from a previous study of the 



magnetization process |jT6[ and ESR-spectrum [[T^ in RbMnBrs. 

The crystal distortions lead primarily to a change in the antiferromagnet coupling be- 
tween chains. But for the purpose of this paper the small distinction between the RbMnBrs 
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magnetic structure and a simple 120°-structure is not so important. Moreover confining 
ourselves by the range of large magnetic fields we assume that the effects of inter-chain ex- 
changes and particularly the distortion of a simple 120° magnetic structure in the absence of 
magnetic fields will give only quantitative corrections not changing the qualitative picture. 
This assumption appears to be true and we estimate these corrections for large magnetic 
fields (see (|l9l-||) below). 



III. EXPERIMENTAL PROCEDURE AND RESULTS 



A vibration-sample magnetometer, similar to the one described in ref. ||T8[, was used 
to measure the magnetization in RbMnBra. A magnetic field up to 6T was generated 
by two superconducting coils. Simultaneous measurements of two mutually perpendicular 
components of the magnetization of the sample, one of which {Miong) is parallel and another 
{Mtr) is perpendicular to the direction of magnetic field were done by using three pairs of 
measuring coils. 

The absolute accuracy of the magnetization measurements was about 5%. The crystal 
was oriented with an accuracy of 1-2°. 

The investigation was performed on single crystals having approximately 1.5 x 1.5 x 
1.5mm'^ volume and approximately 30mg mass in the temperature range 1.7-12K. 

The magnetization component Miong parallel to the field as a function of the field H for 
H\\C and H1.C at T=1.7K is shown in Fig.l. In fields above the critical field {He ^ 3.9T) 
the magnetizations M;^^^ and Mf^^g are not the same, but differ by about 7%. This 
anisotropy of magnetization can not be explained just by taking into account the easy plane 
anisotropy in classical calculations, because it leads to Mi^^g < Ml^^^ while experimentally 
the situation is opposite. Moreover corrections to the classical magnetization due to easy- 
plane anisotropy have to be of order D / J 1% which is much smaller than the observed 
anisotropy (5 — 10%). 

The anisotropy in magnetization is shown on Fig. 2 where the magnetization component 
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Mtr perpendicular to the magnetic field H as a function of the field is plotted at different 
angles between the field and the basal plane and at T=4.2K. Due to the anisotropy of the 
magnetization ai H > He, a non zero transverse-magnetization signal occurs at 7^ 0. 

According to Fig.2 the field dependence Mtr{H) exhibits a small hysteresis. This hys- 
teresis takes place when the projection of the magnetic field on the basal plane reaches the 
value of 3.19 ± O.IOT in increasing fields and 2.67 ± 0.07T in decreasing fields. It corre- 
sponds to the first order phase transition from the incommensurate magnetic structure to 
the commensurate one. 

In Fig. 3 the field dependence of the magnetization component Mtr at different temper- 
atures is shown. At large magnetic fields we can see that Mfr increases with temperature 
supporting the idea that the anisotropy of the magnetization at H > is due to fluctua- 
tions. 

This is even clearer from Fig.4 where the value of Mtr{H = 4.5T, (p = 20°) is shown as 
a function of temperature. 

IV. THEORY 

We start with the microscopic Hamiltonian 

n = + 2J' J2 s^s,+A^ + d ^(5^' -hY.s^ (i) 

i i it 

The intra-chain antiferromagnet exchange constant J is the largest parameter in the Hamil- 
tonian (P and is about two orders of magnitude larger than the exchange constant between 
chains J' and easy-plane single-ion anisotropy D [0, J » J', D and J, J', D > 0. h is the 
magnetic field in the units of energy and it is equal to h = gfisH where if is a magnetic field 
in some standard units, /i^ is a Bohr magneton and g is Lande factor which is equal 2 in the 
compound under consideration. For h = the classical ground state of the Hamiltonian (|I]) 
is the 120° structure with spin vectors forming equilateral triangles in the basal plane. The 
case of a non zero magnetic field was considered classically by Chubukov . For h along the 
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Z (C) axis the transverse components of the spins conserve 120° structure and the classical 
magnetization 

Mn\\c =< >= ^ (2) 

when h < 8JS neglecting by J'/J,D/J (here and there on M means Miong where the 
opposite is not stated explicitly). The case of the magnetic field perpendicular to the Z 
axis is more complicated. For a sufficiently strong easy-plane anisotropy, spins do not leave 
the planes but the 120° structure is not conserved anymore and a spin flip of the two 
magnetic sublattices takes place at the relatively small mag netic field he = VISJJ^. At 
he < h < 8JS we have a coUinear spin structure and magnetization has to be the same (up 
to small y, J- corrections) as in the case of field along Z axis 0. 

In this section we take into account quantum and thermal fiuctuations in the framework 
of linear spin-wave theory and show that the anisotropy of magnetization at h > h^ due 
to these fiuctuations is of the right sign {Mh±c < Mh\\c) and is much larger than the one 
expected from classical calculation. 

In what follows we confine ourselves to magnetic fields h > he where we have a collinear 
spin structure for a magnetic field perpendicular to the C axis of the crystal and we neglect 
the J' term in the Hamiltonian (|l|) in the body of this chapter. This term will give corrections 
to the gap in a spin-wave spectrum of order of J'/ J without changing the qualitative picture 
(see ( pT| - |2^) below). Thus we consider the Hamiltonian of a single magnetic chain instead 
of (|l|) keeping in mind that our results will be applicable only for the h > h^ = v^KLTS^ 
region. We can write Hamiltonians for two directions of magnetic field in the form: 

n^c = 2jy: s^s,+i -hY^st + D Y^{stf (3) 

i i i 

Uh^c = 2JY: S.S.+i -hY^St + D YiStf (4) 

i i i 

Here we slightly changed notations for simplicity. Now magnetic field is always applied 
along Z axis and the easy plane is XY plane for the parallel case and YZ plane for the 
perpendicular one. Let us apply usual Dyson-Maleev transformation 
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S' = 


—S + a^a 


5+ = 


V2Sa+{l 


s- = 


V2Sa 



2S 



(5) 



to (|@) with the axes of quantization z' directed as it is shown in Fig. 5. Choosing the angle 
(j) to cancel linear in a, a"*" terms in Hamiltonians we reproduce the classical results for 
magnetization. The terms of ^ quadratic in a, will give us the linear spin-wave theory 
corrections to these results. We have: 

. , h 



sm I 



sm I 



2S{AJ + D) 
h 



8JS 



(6) 
(7) 



and both Hamiltonians ^ after Fourie transformation can be written in the form 



H = 7 + E 



uji{k)alak - 



with 



DS 



711 = -N{2JS'^(l + 2 sin^ ^) - — ^^^^ ^ + ^-^^ ^in^ , 
ui\\ (k) = 475(1 — sin^ cf) cos k) + DS cos^ 



'^2|| (^) = 475* cos^ cos k + DS cos^ 
7x = -N{2JS\l + 2sm'^(j)) - 



DS 



uJi±{k) = 4J5(1 - sin^ ^ cos k) + DS 
^2±{k) = AJS cos^ (j) cos k — DS 



(8) 



Making Bogolyubov transformation 



ak 



2 Jul - ul 



1/2 / r- 1/2 



we diagonalize (H) and get 



^ = 7 + E 

k 



2JujI-ujI 



1 cJi(A;) 



-k 



(9) 



(10) 



with the spectrum of spin waves given by 



or 



ej(fc) =AJS{l-cosk) [4J5(l + cos20cosfc) + 2L)5cos2 
el{k) = [AJS{l-cosk) + 2DS]AJS{l + cos2(f)cosk) 

Now we can find the magnetization using general relation M = — < 

1 dui 



an 

dh 



> 



M 



5^ r 



2 dh 



or after simple manipulations using (< b^bk > +\] 



M 



dk 


de{ky 




dh 



coth 



1 1 e 
- coth — - 

2 2T 



(11) 
(12) 



(13) 



(14) 



Here the first term describes the classical part of magnetization while the second one comes 
from the contribution of quantum and thermal fiuctuations. 

Let us show now that Eq. (0) leads to the enhanced anisotropy between Mh\\c and 
Mh±c- Consider the case of zero temperature for simplicity. Substituting 1 for coth ^ and 
taking derivatives using (|,^|ll],|l2|) we have 



c 



M, 



hLC 



-S sin I 



— sm ( 



-S sin I 



^ dkl 



[1 — cos A;) 



cosA; + D/AJ 



-n2n2\ 1 + cos 20 cos + Z^/2 Jcos^ 1 + D/AJ 



sm I 



^ dkl ll~cosk + D/2J 



cos k 



(15) 



'-TT 27r 2 V 1 + cos 20 cos 

Now it is clear that only one of two soft modes of spin waves corresponding = vr gives a 
big contribution to the magnetization. This result is quite obvious because it is this mode 
that corresponds to fluctuations of the angle between spins and magnetic field. Another soft 
mode at k = corresponds to asimuthal fluctuations of spins around the direction of magnetic 
fleld and almost does not contribute to magnetization. This mode however contributes very 



8 



much to the value of average spin on the site.Q Expanding the integrands in Eq.([l5D in the 
vicinity oi k = 7T and neglecting D / J m comparison with 1 we get 

dk 1 



M ^ — S" sin + sin j 



27r v/FT^ 



(16) 



with 



9l 



D 



(17) 
(18) 



or taking integrals finally 

Mh\\c 
MhLC 



h_ 
8J 



1 + 



ln( 



+ 



2'kS J 
ln(: 



(19) 
(20) 



2ttS 'I6J2527 

Now let us estimate the corrections to ([T^-^Dp due to 3D (J') effects. These effects give 
corrections of order of J' / J under logarithm in ( |T9| - |20|) and can be calculated in the 
framework of the same linear spin wave theory. Calculations show ||2l| that 3D fluctuations 
will effectively increase the gap of ID system giving instead of ([l9| - pOD : 

h \ 1 / D J'W , , 

The formulas (|2Tl - p2D are good at sufficiently big magnetic fields where fiuctuations and 
effects of three-dimensional ordering are small enough to guarantee the applicability of linear 
spin wave theory for quasi-one-dimensional system. 



^It is seen from this argument that the attempts ||T^ to relate the suppression of magnetization 
with the suppression of average spin by fluctuations are qualitatively wrong because these effects 
arise from the different modes of spin waves. 

^Expressions like In ^ were obtained in |2C] for spin reduction but can not be applied to our case 



(see footnote 1). 



Now from (p^|-p2|) the origin of the enhancement of anisotropy is clear. Big fluctuations of 
the angle between spins and magnetic field give a significant contribution to magnetization 
which depends strongly on the corresponding gap in the spectrum of spin waves due to low 
dimensionality of the system (integral in Eq. (|I6|) would diverge for g = 0). This gap is 
larger for the case of magnetic field along C-axis due to bare easy-plane anisotropy in (|I]). 

V. COMPARISON WITH EXPERIMENT AND DISCUSSION 

To compare our theory with an experiment we need to know the values of the parame- 
ters of the Hamiltonian (|I]) obtained from some independent source. The inelastic neutron 
scattering experiments provide us with such a source. As reported in [|l^ J = 199GHz and 
D = 2.2GHz. These values were found as the best ones to fit the classical formulas for 
spin-wave dispersion to the data of an inelastic neutron scattering. To get the values of J 
and D in we have to take into account the renormalizations by quantum fluctuations. 
Using (4) and (5) from ^ we have: 

J = 186GHz,D = 1.3GHz 

Using these values of J and D and the expressions { ^Vj) and ( p2l) with J' = 0.22GHz from 



He = v48JJ^S^ = 3.9T we have calculated the magnetization of RbMnBrs for both H \\ C 



and H J- C. The range of applicability of ( pij) and (|22|) is approximately H > AT. For 
H < AT the effects of three dimensional ordering and next order corrections to the linear 
spin- wave theory have to be seriously taken into account. 

The magnetizations (|2T| - p2D are plotted in Fig.l together with the experimental data for 
Mh\\c and Mh±c- Also we plotted the magnetizations calculated using classical formulas 
and the values J = 199GHz, D = 2.2GHz. We can see that the formulas (^) and (|2^ ) give 
both the much better absolute values of the magnetizations then the classical expressions 
and describe the anisotropy in magnetization. One can notice also that the experimental 
dependence Mh\\c is a little bit nonlinear with the slope increasing with H (this is not clearly 
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seen from Fig.l but can be checked by fitting the experimental data by linear functions in 
different ranges of H). This can be naturally explained by suppressing fluctuations with 
increasing H and it is seen in (pTI). 

The graphs obtained from (pl|-|2^) in Fig.l lie below experimental points. This probably 
can be corrected by the second order in 1/5 corrections to magnetization which will con- 
tribute with the sign opposite to the sign of the first order corrections. The contribution of 
the higher order terms in l/S" is not very small because the fluctuations are quite big in the 
experimental range of parameters. For example ~ 30% at H = 4.5T. 

VI. CONCLUSION 

The new effect of enhancement of anisotropy in quasi-one-dimensional spin systems due 
to quantum and thermal fluctuations was found. It explains the anisotropy in magnetization 
seen in experiments. The essence of the effect is that the fluctuational part of magnetization 
is sufficiently big and is determined mostly by logarithm of the gap in the spin wave spectrum. 
This gap is anisotropic due to the bare easy-plane anisotropy. The logarithmical dependence 
of the fluctuational part of the magnetization on this gap (essentially the divergence of one- 
dimensional fluctuations) leads to the strong "enhanced" anisotropy of the magnetization. 
The same fluctuations explain a slightly nonlinear character of Mh\\c' dependence and the 
excessive values of magnetization calculated using classical formulas. 

This theory can be applied to the compounds with easy plane anisotropy and quasi-one- 
dimensional magnetic structure such as RbMnBra ,CsMnBr3 and KNiCla 



^For the magnetization data for CsMnBra taken from ||^^ formulas (^l|-22) give even the better 
fit then for RbMnBis in Fig.l. 
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FIGURE CAPTIONS 



Figure 1 The magnetization Miong parallel to the direction of magnetic field as a function of 
field. Circles (o) denote M\\{H) for H \\ C and squares (□) denote M^{H) for H ^C. Solid 
and dashed lines are the results of calculation with parameters J = ISGGHz, D = l.SGHz, 
He = 3.9T, J' = 0.22GHz using formulas (0-g2|) and classical theory § with J = 199GHz, 



D = 2. 2 GHz respectively. 

Figure 2 The magnetization Mfr perpendicular to the magnetic field as a function of H 
with the field making small angles (f with the basal plane, ip = 36° (1), ip = 26° (2), (f = 16° 
(3), ip = 0.5° (4), ^ = -14° (5), V = -24° (6), if = -34° (7). 

Figure 3 The field dependence of the magnetization component Mtj. at different tempera- 
tures with the field making angle v? = 20° with the basal plane. T= 2.4K (1), T= 4.2K (2), 
T= 6.0K (3), T= 7.1K (4). 

Figure 4 The magnetization Mtr{H = 4.5T, (/? = 20°) is shown as a function of temperature. 
The dotted line is a guide to an eye. 

Figure 5 Chosen quantization axes and coordinate systems. 
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